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Abstract 

An approach is presented for the generation of 
two-dimensional, structured, dynamic grids. The grid 
motion may be due to the motion of the boundaries 
of the computational domain or to the adaptation of 
the grid to the transient, physical solution. A time- 
dependent grid is computed through the time integra- 
tion of the grid speeds which are computed from a 
system of grid speed equations. The grid speed equa- 
tions are derived from the time-differentiation of the 
grid equations so as to ensure that the dynamic grid 
maintains the desired qualities of the static grid. The 
grid equations are the Euler-Lagrange equations de- 
rived from a variational statement for the grid. The 
dynamic grid method is demonstrated for a model 
problem involving boundary motion, an inviscid flow 
in a converging-diverging nozzle during startup, and a 
viscous flow over a flat plate with an impinging shock 
wave. It is shown that the approach is more accu- 
rate for transient flows than an approach in which the 
grid speeds are computed using a finite difference with 
respect to time of the grid. However, the approach 
requires significantly more computational effort. 

Introduction 

As computational methods become more efficient 
and computational resources become more powerful, 
the methods of computational fluid dynamics are be- 
ing applied to more complex problems. Among these 
are problems involving the motion of the boundaries 
of the computational domain and problems involving 
transient flow conditions with solution adaptive grids. 
These generally require a dynamic grid. Much work 
has been done on methods for the generation of static 
grids with proper qualities. A dynamic grid should 
maintain these same qualities; however, this topic has 
not received much attention. 

We present an approach for the generation of dy- 
namic grids for which the qualities of grid smoothness, 
orthogonality, and cell volume adaptation are main- 
tained. The dynamic grids are generated through the 
time-integration of the grid speeds computed from a 


system of grid speed equations. The grid speed equa- 
tions are derived from the time- differentiation of the 
grid equations. This approach was first presented by 
Hindman et al. [1] for an application involving the mo- 
tion of a shock-fitted boundary. Holcomb and Hind- 
man [2] extended the approach for the adaptation of 
the grid to the physical solution. The previous cited 
work involved grid equations based on the Poisson 
equations and did not account for grid orthogonality. 
The present approach uses a system of grid equations 
based on the variational approach of reference [3]. 

An alternative approach for computing the grid 
speeds is to use a time-difference of the grid. This ap- 
proach is examined. It will be shown that solving the 
grid speed equations reduces the grid-motion-induced 
errors for the inviscid, unsteady flow in a nozzle. 

The objectives of this paper are to present the 
system of grid and grid speed equations for two- 
dimensional space, present an explicit dynamically 
adaptive grid method, and demonstrate the behav- 
ior of the method through applications. These ap- 
plications include a simple problem involving a mov- 
ing boundary, an inviscid flow through a converging- 
diverging nozzle during startup, and a viscous flow 
over a flat plate with an impinging oblique shock. The 
last two flows involve flow conditions which are ini- 
tially transient; however, after time, the flow condi- 
tions and the flow reach a steady-state. 


The Grid Equations 

The grid equations follow the variational approach 
of Brackbill and Saltzman [3] which derive the grid 
equations from the Euler-Lagrange equations, 


(d_ _d_ d d d \ 

V <9 f dx^ + dr) dx v dx ) ~ ^ 


(1) 


and 

(d_ _d_ d d 8 \ T 

\<9£ + dr / dy n dy) ~ °‘ ^ 

The LagrangianX is defined so as to include measures 
of smoothness, orthogonality, and cell volume adapta- 
tion and is of the form 


L — Ls 4- Lq + La. (3) 
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The expressions for Lo and L A presented here are 
a modification of those presented in reference [3] to 
nondimensionalize the Lagrangians. The Lagrangians 
are of the form 


= A rf r n 


(4) 


L o = Ao (±A)\ (5) 

and 

L a = X A (W J f / K\ (6) 

The Xs, A o, and A^ are specified constants which 
weigh the importance of each Lagrangian. The cell vol- 
ume adaptation Lagrangian L A is essentially a state- 
ment of the equidistribution of the product WJ over 
the grid. The W is a weighting function usually di- 
rectly related to the gradients in the flow solution. The 
J is the Jacobian and is a measure of the cell volume. 
The K is the average value of WJ over the grid. 

The system of grid equations are derived by in- 
serting the expressions for the Lagrangians into the 
Euler-Lagrange equations and is expressed in the form 


G( f) = Artf + Br^ + C r^ + Drz + Err, = 0 (7) 

where r T = [x, y ] . The coefficient matrices are sums 
of the form 

A = As + Ao -f- A a ■ 

The B, C, D , and E matrices are of the same form. 
The expressions for the smoothness matrices are 


A s = ^[S], (8) 

Bs = ~j X 3 S ~ [5] , (9) 

c s = Ay t s] , ao) 

and Ds and E s are the null matrix [0]. Also, 

M-eVj' 

The remaining variables are defined as 


J = Vti - Xr, Vt, ( 12 ) 

<X = *? + !/?. (13) 

0 = X£ x v + yt y ni (14) 

7 = x 2 + y 2 , (15) 

= y 2 + y 2 , (16) 

€ X£ y^ -4- x 1 7 t/fj , (17) 


and 


p = x\+ *5. 


(18) 


Ao = 


The expressions for the orthogonality matrices are 
A ooi 


J 4 


2 Pihyt + ttyj ip 

fp 20x n x^ + ax | J 


i (19) 


\ o {20 2 + ay) 

° = -p — 


2 y^y v - <P 

- <p 2x^x r) J 


( 20 ) 


Co = 


Xpy 

J A 


7 y 2 + 2/%^ 


yx 2 + 20XfXt, J 


( 21 ) 


and Do and Eo are the null matrix [0]. The <j>, ip, and 
k are defined as 


<P = a* y n + x n y f} (22) 

xP ~ - ( y Xr} y v + 2 a x$ y$ ), (23) 

and 

k = - ( a vs + 2 7 x n y v ). (24) 

The expressions for the cell volume adaptation 
matrices are 


A a = 


Aa ^ 2 

K 2 


y 2 v 

Xr/yri 


Xt] yr) 
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B a = 


A A W 2 
K 2 


- 2 y^y n <p 
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C A 
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D a 


A A JWW n 

K 2 


0 - 1 
1 0 


and 


E _ A A JWW i 

K 2 


0 1 

- 1 0 


(25) 

(26) 

(27) 

(28) 

(29) 


The Grid Speed Equations 

The grid speed equations are obtained from the 
time-differentiation of the grid equations. Applying 
the chain rule of differentiation and collecting terms 
results in the grid speed equations of the form 


C r (x, y) — A* + B* -f C* z vv 

+ D* Z{ + E * z n + f* = 0 (30) 
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where the grid speed vector is defined as 
z — [ x T , 2/7- ] • 

This involves performing the time differentiations A r , 
B T) C t , D t , and E r and factoring out z. The proce- 
dure is straight-forward; however, the algebra is quite 
involved. The details are presented in reference [4]. 
One observation is that the time-differentiation of the 
matrices only result in first-order time derivatives of 
the grid and so 

[A*] = [A], [B*] = [£], and [C*j = [C]. 

The Weighting Funct ion 

The weighting function W is defined in the form 

W = A 0 + Ai | V/(f/)| + A 2 exp ( A 3 |V/(U)| ) . (31) 

The V f(U) is the gradient with respect to the compu- 
tational coordinates £ and t) of a function /, such as 
density, of the physical solution. The function / may 
be discontinuous, thus, an elliptic smoothing is used 
to make / a continuous function. The Ao, Ai, A 2 , and 
A3 are specified constants which influence the spatial 
character of W . 

The grid equations require the derivatives and 
Wrj. The grid speed equations also require the deriva- 
tives W T , Wf T , and W VT . These are computed by ap- 
plying the chain rule and using second-order spatial 
and first-order temporal finite-differences. 


applied and an expression can be developed for the 
change in the curve-length coordinate s needed for or- 
thogonality for the iteration. For the j — 1 boundary, 
this expression is 


(A*k? = 


(^)»,i ' [ 4 ^~, 2 ~ 3f-,i - r it3 ] 
3 (r* -r s ) iA 


, (33) 


where m is the iteration index. Thus 


s m+1 = + (As) m (34) 

and the new boundary point grid point is 

r B m+1 = fflt ( * m+1 ) (35) 

where the subscript B denotes a boundary point and 
Ffp. is the curve fit of the boundary. 

The grid speeds at the boundaries are of the mixed 
type to impose a prescribed boundary motion and to 
impose the orthogonality condition at the boundary. 
The boundary condition for imposing orthogonality is 
derived from the time differentiation of the grid or- 
thogonality condition 

Q 

g;(r r r-) = 0. (36) 

Using the parameterization of the boundary and a lo- 
cal linearization, an expression for the the time deriva- 
tive of the curve-length coordinate can be obtained. 
For the j = 1 boundary, the expression is 


Solving the Grid and 
Grid Speed Equations 

The grid equations are non-linear partial differen- 
tial equations while the grid speed equations are linear 
partial differential equations. This is attractive since 
linear equations are easier to solve than non-linear 
equations. 

The smoothness terms in the grid equations form 
a system of elliptic equations. The orthogonality and 
cell volume adaptation terms form a system of equa- 
tions which is not elliptic and can produce non-unique 
solutions. However, the combined system of equations 
remains fairly elliptic and iterative methods for solving 
elliptic equations can be applied with success. 

The grid and grid speed equations are discretized 
using second-order, central differences and solved us- 
ing the Gauss-Seidel point relaxation method. 

Neumann boundary conditions are applied to im- 
pose orthogonality of the grid lines at the boundaries. 
The orthogonality condition is 

^ • r v = 0. (32) 

The boundary is parameterized according to the curve- 
length along the boundary. A local linearization is 


«5r = 


' ( 4 n, 2 ~ 3rj,i - f, |3 ) + r s • (4 z i<2 - z it 3 ) 
3(r s • r s ) 


Thus for each iteration, a new s T is computed and 


(37) 


s m+l _ s m + Ar ^ 

and from the curve fit of the boundary, 

?B m+1 = % ( s m+ ‘ ) • (39) 

The mixed boundary condition is imposed as 

z b{^) = z B(motion)(t) + ( f g ~ r g 1 ) / A T, (40) 

Solving the grid speed equations for the grid 
speeds and then integrating to obtain the grid is ex- 
pected to generate a grid which satisfies the grid equa- 
tions. However, small deviations may develop. One 
way to ensure that the grid continues to satisfy the 
grid equations is to rewrite the grid speed equations 
in the form of a first-order homogeneous dynamic sys- 
tem [2] as 

G t + Ac G = 0, (41) 

where A c is a prescribed damping constant which 
damps deviations from the grid equations. 
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Coupling of the Flow 
and Grid Equations 

The flow equations and the grid speeds are in- 
tegrated in time using an explicit, two-stage Lax- 
WendrofF method [5] of the form 

= *.-3 + (42) 

+Ar (43) 

and 

= \ { + *ij* } ■ (44) 

For the grid, $ = r and ■= £. For the flow equations, 
4* — U and — R. The U and R denote the alge- 
braic vectors of the generalized conservative variables 
and fluxes defined from a cell- vertex, finite-volume ap- 
proximation of the unsteady Navier-Stokes equations. 
Thus, 

u = U V, (45) 

where V is the area of the planar finite-volume cell 
and £/ is the algebraic vector of conservative variables. 
The R is the flux residual defined as 

Rij — ^»'+l/2,j — Fi-l/2,j + Fi,j + l/2 - -FjJ-1/2) (46) 

where the flux vector for a cell face is defined as 

F = H - hdS. (47) 

The H is the flux dyadic, which for a time-varying 
control volume is 


H - F-zU. (48) 

The F is the Cartesian flux dyadic for the two- 
dimensional, unsteady Navier-Stokes equations. The 
flow equations are complete with Sutherland’s formula, 
the definition of the Prandtl number, and the assump- 
tions of a perfect gas (air) and laminar viscous flow. 

The time step is computed using the CFL condi- 
tion. The inviscid fluxes are computed using the Roe 
flux-difference splitting with a TVD limiter as pre- 
sented in reference [5]. The viscous fluxes are com- 
puted using differences and averages computed at the 
cell faces. Characteristic boundary conditions are used 
to compute the solution points on the boundaries. 

The V needed to decode U is computed from the 
grid conservation law as discussed in references [6] and 
[7]. The grid conservation law relates the change in 
volume of the cell to the motion of the cell faces and 
is derived from the flow integration equations with the 
assumption of a uniform solution for the conservative 
variables. The grid conservation law follows the form 
of equations (42) to (44) with $ = V and ^ = Zwhere 
Z is the vector sum of speeds of the cell faces. 

The flow and grid equations are coupled through 
the grid speed terms in the flow equations and the W 


function in the grid speed equations. The two-stage 
approach advances the solution and grid in time and 
space in ideally second-order accuracy. The steps in 
the explicit, dynamically adaptive grid method using 
the grid speed equations include: 

1) Establish an initial grid F 1 satisfying the grid equa- 
tions and a consistent initial solution U n . 

2) Solve the grid speed equations for F 1 . 

3) Compute or specify the time step, At. 

4) Compute U* for the first stage, 

(IF = U$ - A tR$. (49) 

5) Solve the grid conservation law, 

ViJ = Vtf-ATZJ. (50) 

6) Compute the first-stage flow solution, 

u h = Vtj / ViJ. (51) 

7) Compute the flow boundary conditions. 

8) Integrate the grid speeds to obtain the grid, 


= r. 




+ At zF- . 


(52) 


9) Recompute the cell volumes and cell-face areas. 

10) Solve the grid speed equations for £ *. 

11) Compute U ** for the second stage, 


U** = U*- 


At R. 


»..7 ' 


12) Compute 


frn+l F 

tj - 2 


u?j + t//; 

13) Solve the grid conservation law, 


(53) 


(54) 


Vi" = ViJ ~ A ri { ]. (55) 

14) Compute the volume at the n + 1 time level, 

V <J +1 = 5 ( Vij + V,r ) . (56) 

15) Compute the flow solution, 

KF = / v ij +i . (w) 

16) Compute the flow boundary conditions. 

17) Integrate the grid speeds to obtain the grid, 


f n+l 




hj 


+ 


1 

2 


Ar 


( £ n 




+ ?<j ) 


(58) 


18) Recompute the cell volumes and cell-face areas. 

19) Solve the grid speed equations for £ n+1 . 

20) Go to step 3 and repeat until the final time. 
Computing the grid speeds using a time-difference 

of the grid requires solving the grid equations at each 
stage prior to the differencing. 
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Figure 1: The grid and grid speed vectors for the ro- 
tation of the lower boundary of a unit square 

Results 

A Simple Dynamic Boundary 

The dynamically adaptive grid method was ap- 
plied to a uniform grid on a unit square domain in 
which the lower boundary was prescribed to rotate at 
a rate of 20 degrees/second about the lower left cor- 
ner. This problem tested the ability of the method to 
compute the grid speeds and integrate them in time 
to create a grid which continuously satisfied the grid 
equations, and so, maintain the desired grid qualities. 

The time integrations were performed for 50 time 
steps with a time step of A r = 0.02 seconds to a final 
time of 1 second. The bottom boundary rotated 20 de- 
grees about the lower left corner. The grid parameters 
were (As,Ac>, A^) = (1.0, 1.0, 1.0) and A c = 30.0. 

For an iteration tolerance of 1,0® 10~ 8 , the com- 
putation required 17308 iterations of the grid speed 
equations. The average error was 0.0396% and the 
maximum error was 0.8176%. The errors are the dif- 
ferences in grid point locations between the final grid 
obtained from the time integration and the grid which 
satisfied the grid equations at the final time. The er- 
rors are expressed as percentages of the length of the 
uniform grid spacing. Figure 1 shows the grid and 
grid speed vectors at the final time. The final grid 
shows that the grid maintained the desired grid quali- 
ties, specifically, the orthogonality of the grid lines to 
the lower boundary. 

The effectiveness of the grid control law damping 
factor, Ac, is shown in figure 2. The average error 
was reduced for increased values of Ac. The fact that 
the curve bottomed out indicates that one should be 
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Figure 2: The effectiveness of the Ac damping factor 

able to pick a value of Ac and obtain a good reduc- 
tion in error without having to search for an optimum 
value. However, while the errors decreased, the equa- 
tions became stiffer, and more iterations of the grid 
speed equations were required. The effectiveness of 
Ac remained consistent with other test problems, and 
so, these results guided the selection of Ac for the flow 
problems. 

The pointwise iteration of the grid speed equa- 
tions demonstrated excellent convergence characteris- 
tics. When the iteration tolerance was increased to 
1.0® 10 -4 , only 304 iterations were required and the 
average error and maximum error only increased to 
0.1065% and 2.8635%, respectively. When the num- 
ber of iterations per call was limited to 5 iterations, 
the total number of iterations was 332 and the average 
and maximum errors were only 0.1023% and 2.0723%, 
respectively. This provided some guidance for setting 
the number of iterations and tolerances appropriate 
for the efficient solution of the grid speed equations. 

Inviscid Startup Flow in a Nozzle 

The dynamically adaptive grid method was ap- 
plied to the computation of the inviscid startup flow in 
a converging-diverging (CD) nozzle. Figure 3(a) shows 
the shape of the nozzle. The flow was initially station- 
ary at a total pressure of p t = 104074.6 N/m 2 and a 
total temperature of T) = 252.9 K . The exit pressure 
was then decreased from the total pressure value to a 
value of p ex i t = 84000 N/m? over a time interval of 
10% of a reference time. The reference time was de- 
fined by the ratio L/Vi n f\ ow where L is the length of 
the nozzle ( L = 10 meters) and Vi n fi ow is the inflow 
velocity (Vi n fi ow = 75.92 meters/sec). A two-point, 
cubic spline variation with zero end-slope conditions 
was used for the time variation of the exit pressure. 
After the transition time interval, the exit pressure 
was held fixed. The flow was computed until a final 
time of 0.3 seconds. All computations were performed 
on the NASA Lewis Cray XMP. 

The grids consisted of 95 grid points in the stream- 
wise direction and 15 grid points in the traverse di- 
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Figure 3: The grid and density contours for the static 
grid computation for the CD nozzle 

rection. The static grid, shown in figure 3(a), was 
computed using the grid equations with parameters 
= (1.0, 1.0, 0.0). The dynamic grids 
were computed with grid parameters (A 5 ,Ao,A, 4 ) = 
(1.0, 1.0, 0.5). The function / used in computing the 
weighting function W was the density. The grid adap- 
tation parameters were A 0 = 1.0 and A x = 500. Thirty 
smoothing passes were used to smooth the density so- 
lution prior to computing the weighting coefficient W 
and its derivatives. 

The flow solution on the static grid was computed 
using a CFL number of u — 0.7 and required 3016 time 
steps and 111 CPU seconds to converge to a residual of 
1.308#10“ 3 after a maximum of 6.284#10 -1 . Figure 3 
shows the grid and density contours at the final time. 

The flow solution on the dynamically adaptive grid 
using time-differenced grid speeds was computed us- 
ing a CFL number of u = 0.6 and required 4423 time 
steps and 563 CPU seconds to converge to a resid- 
ual of 1.198# 10 -3 . The grid equations were iterated 
with a tolerance of 1.0# 10 -4 with a maximum of 3 
iterations allowed for each stage of the Lax-Wendroff 
method. One would expect that integrating the grid 
speeds prior to solving the grid equations would pro- 
vide a good first guess of the new grid. Experience has 
shown that for this problem, performing the integra- 
tions caused the boundary of the grid to deviate sig- 
nificantly from its true geometry and introduced per- 
turbations into the flowfield. 

An instability occurred when applying the invis- 
cid wall boundary conditions for a dynamic grid. The 
grid points on the boundary are constrained to move 
along the boundary. Since the boundary conditions 
act only on information normal to the boundary, the 
grid speeds do not influence the boundary conditions. 
However, the numerical implementation of the bound- 
ary conditions may result in some influence from the 
grid point motion. For the dynamic grid computa- 
tions, the grid speeds were neglected in the numerical 
boundary condition computations. 

The flow solution on the dynamically adaptive grid 


using the grid speeds computed from the grid speed 
equations was computed using a CFL number of v — 
0.6 and required 4520 time steps and 944 CPU sec- 
onds to converge to a residual of 9.633#10“ 4 . The 
grid control law damping factor was set to a value of 
A c — 50.0. The grid speed equations were iterated 
with a tolerance of 1.0#10“ 4 with a maximum of 3 
iterations allowed for each stage of the Lax-Wendroff 
method. At the start of the computation, the grid 
speed law required all 3 iterations. As the steady-state 
flow solution was reached, the residual fell below the 
tolerance and less iterations were required until only 
1 iteration was required per stage. Figure 4 shows the 
grids and density contours at various times. As can be 
seen, the resolution of the shock is clearly improved 
over the computation on the static grid. 

The variation of the static pressure in the stream- 
wise direction for all of the computed solutions com- 
pared very well with quasi-one-dimensional, inviscid 
theory. The dynamically adaptive grid solutions show 
increased resolution of the shock. The TVD limiter 
eliminated any oscillations near the shock. 

The convergence histories were similar for each 
computation. The convergence history for the dy- 
namic grid computation which used the grid speed 
equations showed some occasional sharp spikes which 
may indicate movement of a grid point across the 
shock. 

Using a dynamically adaptive grid required signif- 
icantly more computational effort than a static grid. 
The dynamic grid computation using time-differenced 
grid speeds required 3.46 times more CPU time per 
step than the static grid computation, while the dy- 
namic grid computation using the grid speed equa- 
tions required 3.07 times more CPU time. Using the 
grid speed equations was more efficient for this prob- 
lem because linear equations are solved rather than 
non-linear grid equations. Using a time difference re- 
quired 3% more memory than the static grid compu- 
tation, while using the grid speed equations required 
20% more memory. 

The existence of grid-motion-induced errors was 
examined by comparing the time variation of the static 
pressure at the location # = 8.0 meters in the nozzle 
for the static grid computation with the dynamic grid 
computations. As shown in figure 5, there is less error 
when the grid speed equations are used than when a 
time difference is used to compute the grid speeds. 

Shock /Boundary- Layer Interaction 

The dynamically adaptive grid method was ap- 
plied to the computation of a viscous flow in which an 
oblique shock impinges on a laminar boundary layer 
along a fiat plate at # = 1.0 meters and causes a small 
flow separation region. This problem was investigated 
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Figure 4: The grids and density contours at various 
times for the dynamically adaptive grid 



Figure 5: The time history of the pressure at the lo- 
cation x = 8.0 for the startup of the CD nozzle 

experimentally by Hakkinen et al. [8] and is used ex- 
tensively to validate CFD methods. The freestream 
conditions prior to the shock are a Mach number of 
2.0, a Reynolds number of 2.96aT0 5 based on a length 
of 1.0 meters, and a shock angle of 32.585 degrees. 

The flow domain extends from x = -0.2 to x — 2.0 
meters and from y = 0 to y — 1 meters. The grids 
for the flow computations followed the suggestions of 
Liou [9] based on triple-deck theory. The grids con- 
sisted of 75 streamwise and 65 transverse grid points. 
The initial solutions included the incident and reflected 
oblique shocks whose properties were computed from 
shock theory. The flow solutions were computed with- 
out TVD limiters. The velocity at the plate was de- 
creased during the flow computation from the inviscid 
flow velocity to zero velocity over a time interval of 
3.95a? 10 — 3 seconds. 

The static grid was generated using an interactive 
grid generator with a normal grid spacing at the plate 
of 1 .6a? 10“ 4 meters and clustering to increase resolu- 
tion along the inflow boundary and plate. The grid 
is shown in figure 6. The flow solution on the static 
grid was computed using a CFL number of v - 0.8 
and required 90000 time steps and 11519 CPU seconds 
to converge to a residual of 9.8xl0“ 5 after reaching a 
maximum of 5.0. Figure 6 shows the pressure contours 
for the steady-state solution. 

The dynamic grids were computed with grid pa- 
rameters (As, Ao, A^) = (1.0, 1.0, 1.0). The function / 
used in computing the weighting function W was the 
Mach number. The grid adaptation parameters were 
Ao = 1.0 and Ai == 500. Thirty smoothing passes were 
used to smooth the Mach number prior to computing 
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W and its derivatives. The initial grids were obtained 
by solving the grid equations for the given adaptation 
parameters with the flow solution being defined by in- 
viscid, oblique shock theory. Since no boundary layer 
present, there was no clustering normal to the plate; 
the clustering developed as the boundary layer devel- 
oped. Figure 7(a) shows the initial grid. 

The flow solution on the dynamically adaptive grid 
using the time-differenced grid speeds was computed 
using a CFL number of u — 0.7 and required 6658 time 
steps and 2118 CPU seconds to converge to a residual 
of l.OxlO -4 . The CPU time per step averaged 2.49 
times the time required for the static grid. 

The flow solution on the dynamically adaptive grid 
using the grid speed equations used a CFL number of 
v = 0.7 and required 8392 time steps and 2703 CPU 
seconds to converge to a residual of l.OxlO -4 . The grid 
control law damping factor was set to a value of Ac = 
50.0* The iterations of the grid speed equations were 
limited to 10 iterations per stage and the grid residual 
tolerance was set to 5.0xl0 -4 . At the start of the 
computation, all 10 iterations were required. As the 
steady-state flow solution was reached, the residual fell 
below the tolerance and less iterations were required 
until only 1 iteration was required per stage. The CPU 
time per step averaged 2.52 times the time required 
for the static grid. Figure 7(b) shows the adapted grid 
at the final time and figure 7(c) shows the pressure 
contours. 

The adapted grid seems to have improved the res- 
olution of the shocks. Figure 8 shows the comparison 
of the pressure coefficients along the plate. Figure 9 
shows the comparison of the skin friction coefficient 
along the plate. The static grid results compare well 
to the experimental data and are consistent with nu- 
merical results of other researchers. The dynamic grid 
predicted a smaller separation region. This is due to 
the improper resolution of the boundary layer. The 
spacing of the dynamic grid next to the plate was 
approximately 1.6xl0 -3 meters. At x — 1.0 meters, 
this corresponds to [u/u^ = 0.038, compared to 

|u/uoo| = 0.0025 for the static grid. The Uoo is the 
velocity just outside the boundary layer. 

Obtaining a smaller grid spacing normal to the 
plate for the dynamically adaptive grid method relies 
on emphasizing adaptation over smoothness. This in- 
volves increasing the values of A,t , Aj , A 2 , and A3. At- 
tempts to reduce the spacing have yet been unsuccess- 
ful and is under further study. 

Conclusions 

An explicit, dynamically adaptive grid method has 
been presented and demonstrated for two-dimensional 
problems. The dynamic grids maintained the desired 
grid qualities while allowing grid motion due to mo- 
tion of the boundaries and to adaptation of the grid 



Figure 6: The grid and pressure contours for the 
shock/boundary-layer interaction computation on the 
static grid 


to solution gradients. The results suggest that using 
the grid speed equations to compute the grid speeds 
reduces the grid-motion-induced errors more than us- 
ing a backwards time-difference of the grid. The grid 
adaptation improved the shock resolution; however, 
the resolution of boundary layers needs to be im- 
proved. Using a dynamically adaptive grid required 
significantly more computational effort than using a 
static grid. The problems presented attempted to 
demonstrate the approach, but static grids are prob- 
ably adequate for these problems. Future efforts will 
focus on applying the method to truly unsteady flow 
problems and problems involving boundary motion. 
The method is applicable to three-dimensional prob- 
lems and development has begun. 
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